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^ ! Abstract 
u'. . 

D . We investigate the role that baryon number violating interactions may play in B phe- 

Q I nomenology. Present in various grand unified theories, supersymmetric theories with i?-parity 

^ ■ violation and composite models, a diquark state could be quite light. Using the data on B 

CN . decays as well as B — B mixing, we find strong constraints on the couplings that such a light 

I ! diquark state may have with the Standard Model quarks. 
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O ; 1 Introduction 

^ : 

Q-ii It is now widely accepted that the Standard Model (SM), despite its great success, is only an 
Qh' effective theory. The ills plaguing it may be cured only in the context of a more fundamental 
^ ■ theory operative at higher energies. The quest to find such a theory has, over the years, inspired 
• • , many a model going beyond the SM. Two of the most attractive classes of such models are 
\ those incorporating grand unification Q and/or supersymmetry [||. The exact nature of such 
^ ■ a theory, however, is a matter of intense debate, occasioned, not in the least, by the absence 
of any experimental signature yet. It is not surprising, thus, that the search for new physics 
effects constitutes a major component of research in high energy physics. Such efforts can be 
broadly classified into two categories. On the one hand there are the direct searches typified by 
high energy collider experiments where new particles are sought to be produced on-shell and 
detected through their subsequent decays. The other approach concentrates on indirect effects 
as can be deduced from possible deviations from the SM predictions for low-energy observables. 
In this article we shall focus on one such set of low and intermediate energy experiments. 

The next decade will see a blossoming of experimental facilities planning to explore B — B 
mixing as well as 5-meson decays with greater accuracy and for an increasing number of 
different final states. In light of these upcoming experiments (CLEO, BaBar, BELLE, HERAB, 
BTEV and LHCB), it is of importance to examine their sensitivity to new physics beyond the 
SM. 

In this paper, we investigate the possible infiuence that a baryon number violating inter- 
action may have on B- phenomenology. Within the SM, baryon {B) and lepton (L) number 
conservations come about due to accidental symmetries. In other words, such conservations are 
not guaranteed by any principle, but are rather the consequences of the choice of the particle 
content 0. In extensions of the SM, such an accidental occurrence is obviously not guaranteed. 

^Indeed, nonperturbative effects within the SM itself do break B + L symmetry. 
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For example, even in the simplest grand unified theories (GUTs), both the gauge and the 
scalar sector interactions violate each of 13 and L. The corresponding particles, namely the 
diquarks Q and leptoquarks have been studied in the literature to a considerable extent. 

Simultaneous breaking of both B and L symmetry is obviously a recipe for disaster as this 
combination is more than likely to lead to rapid proton decay. Within GUTs, gauge boson- 
mediated proton decay could be naturally suppressed by postulating the symmetry breaking 
scale to be very large. However, there do exist a class of GUTs |5|, where the next set 
of thresholds need not be very high and i?-violating gauge particles can be relatively light. 
Proton decay, however, remains suppressed on account additional symmetries in the theory. 
Suppression of the scalar mediated contribution to proton decay in a generic GUT, on the 
other hand, is easier to obtain: the particle content can be so chosen that there is no diquark- 
leptoquark mixing, at least as far as the light sector is concerned. 

In the case of the Minimal Super symmetric Standard Model (MSSM), though, we do not 
have the option of demanding the 'offending' fields (the supersymmetric partners of the SM 
fermions) to be superheavy. Ruling out the undesirable terms necessitates the introduction of 
a discrete symmetry, R = (— l)3(S-L)+2S (^-^^^i^j^ § denoting the spin of the field) Q. Apart 
from ruling out both B and L violating terms in the superpotential, this symmetry has the 
additional consequence of rendering the lightest supersymmetric partner absolutely stable. 
However, such a symmetry is ad hoc. Hence, it is of interest to consider possible violations 
of this symmetry especially since it has rather important experimental consequences, not the 
least of which concerns the detection of the supersymmetric partners. 

It can thus be argued that, in such models as well as in models of compositeness Q, it is 
quite likely that baryon number violating interactions may not be suppressed too severely. Even 
more interestingly, such processes may be mediated by relatively low-lying states, generically 
called diquarks. 

This paper is organised as follows: Section 2 constitutes a brief review on diquarks. Section 
3 deals with hadronic B decays. In section 4, we concentrate on B — B mixing. Section 5 
contains the numerical results. We conclude in section 6 with a summary and outlook. 



2 Diquarks: a brief review 

In this section we shall briefly examine all possible tree-level i?-violating couplings involving 
the SM quarks. We shall adopt a purely phenomenological standpoint without any particular 
reference or prejudice to the origin of such couplings or states. A generic diquark is a scalar or 
vector particle that couples to a quark current with a net baryon number B = ±2/3. Clearly, 
under SU{3)c, it may transform as either a triplet or a sextet. For scalars, the Yukawa term 
in the Lagrangian can be expressed as 

^SD = h\fq^PL,Rqj<^A + h.c, (la) 

where i, j denote quark flavours, A denotes the diquark type and Pl,r reflect the quark chirality. 
Standard Model gauge invariance demands that a scalar diquark transforms either as a triplet 
or as a singlet under SU{2)l. For a vector diquark, on the other hand, the relevant term in 
the Lagrangian can be parametrized as 

CvD = ^\fq'lt.PL,RqjV^ + h.c. (lb) 

with Va transforming as a SU{2)l doublet. The full list of quantum numbers, for either case, 
is presented in Table |l[ Clearly, the couplings h'lj\ h[^^ and /ijj^ must be symmetric 
under the exchange of i and j while h^^\ h^^\ h[^^ and h^^^ must be antisymmetric. For the 
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other couplings, viz. , , and ^Ij , no such symmetry property exists, hereafter, we 
assume these coupHngs to be real 0. Note that the quantum numbers of $2,4,6 as well as those 
of V24^ allow them to couple to a leptoquark {i.e. a quark-lepton) current as well. Clearly, the 
non-observance of proton decay implies that such L-violating couplings must be suppressed 
severely. 

It should be noted that we are not demanding that the vector diquarks correspond to 
some gauge theory. While it might be rightly argued that a theory with non-gauged vector 
particles is non-renormalizable, one should keep in mind that such states may well be there in 
an effective theory. Since we would be studying the phenomenological implications only at the 
lowest order of perturbation theory, renormalizability is not an issue here. 



Diquark Type 


Coupling 


SU{2,)c X SU{2)l X U{1)y 


$1 




(6, 3, -|) 


$2 


h'^\QLirQLj<^2 


(3, 3, -|) 


$3 


[t^hQLiTQLj + hf^HuRiYdRj] 


(6, 1, -I) 


$4 


[h^\QLirQLj + h^\uRirdRj] CD4 


(3, 1, -1) 


$5 




(6, 1, -1) 


$6 


h\f {uRiYuRj<^G 


(3, 1, -I) 


$7 


h^\dRi)^dRj^7 


(6, 1, 1) 


$8 


hf{dRi)^dRj^S 


(3, 1, 1) 


Vif" 


4fiQuri^.dR,V,^ 


(6, 2, i) 




^f^iQuTl^dR.V^ 


(3, 2, i) 






(6, 2, -|) 




^f{QL^rl^.uR,Vt 


(3, 2, -f ) 



Table 1: Gauge Quantum Numbers and Yukawa Couplings of Diquarks {Qem = 73 + ^)- 

We now turn to the MSSM, where both B- and L-violating terms are allowed, in general, 
by supersymmetry as well as gauge invariance. As stated earlier, catastrophic rates for proton 
decay can be avoided by imposing a global Z2 symmetry under which the quark and lepton 
super fields change by a sign, while the Higgs super fields remain invariant. However, since such 
a symmetry is entirely ad hoc within the purview of the MSSM, it is conceivable that this 
i?-parity may be broken while keeping either of i? or L intact. In our study, we shall restrict 
ourselves to the case where only the i?-violating terms are non-zero. Such scenarios can be 
motivated from a class of supersymmetric GUTs as well [^. The corresponding terms in the 

^The extension to complex couplings is straightforward. The imaginary parts, however, can be better con- 
strained from an analysis of the CP violating decay modes. 
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superpotential can be parametrized as 



= Ar,.,f7)j^ijDi (2) 

where TJ\^ and Z))^ denote the right-handed up-quark and down-quark superfields respectively. 
The couphngs A^^-^ are antisymmetric under the exchange of the last two indices. The corre- 
sponding Lagrangian can then be written in terms of the component fields as: 

= A-.-fc (<44 + + u\d]d'C) + h.c. (3) 

Thus, a single term in the superpotential corresponds to two of type li'^^ and one of type lif^ 
diquark interactions. 

The best direct bound on diquark type couplings is derived from an analysis of dijet events 
at the Tevatron Considering the process qiq^ $a QiQj, an exclusion curve in the 
{m^j^,h[^^) plane can be obtained from this data. A similar statement holds for the vector 
particles as well. Two points need to be noted though. At a pp collider like the Tevatron, the 
uu and dd fluxes are small and hence the bounds are relatively weak. This is even more true 
for quarks of the second or third generation (which are relevant for the couplings that we are 
interested in). Secondly, such an analysis needs to make assumptions regarding the branching 
fraction of V^a) into quark pairs, a point that is of particular importance in the context of 
-R-parity violating supersymmetric models. 

There also exist some constraints derived from low energy processes. Third generation 
couplings, for example, can be constrained from the precision electroweak data at LEP or, 
to an extent, by demanding perturbative unitarity to a high scale Couplings involving the 
first two generations, on the other hand, are constrained PI by the non-observance of neutron- 



antineutron oscillations or from an analysis of rare nucleon and meson decays |11, 1^. While 



many of these individual bounds are weak, certain of their products are much more severely 



constrained by the data on neutral meson mixing and CP-violation in the i^T-sector |13|. It 
is our aim, in this article, to derive analogous but stronger bounds. 

At energy scales well below the mass of the diquark, the latter can be integrated out and 
effective four quark operators obtained. In Table ^, we list these for each diquark type. A few 
points should be noted here: 

• we have not displayed the operators resulting from and as these do not contribute 
(at the lowest order) to either B decays or B — B mixing; 

• for convenience's sake, we have Fierz-rearranged the operators and, in the process, ex- 
changed the charge-conjugated fermion fields (which come in naturally) for their non- 
conjugated counterparts; 

• within a diquark multiplet, we have assumed all the fields to be mass degenerate since 
large splittings within a multiplet are anyway disfavoured by LEP data; 

• We have neglected the evolution of the diquark mediated effective four-quark interac- 
tions from the electroweak scale down to B meson scale through renormalisation group 
equations; 

• we have not displayed the extra color factors that appear on account of the diquark 
states being coloured objects. The said factors can be determined by reexpressing four 
quark operators of the forms (3c 3c) i and (6c ® 6c) i in terms of the corresponding 



^Although many of these analyses have been done for the case of i?-parity violating models, clearly similar 
bounds would also apply to nonsupersymmetric diquark couplings as well. 
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Diquark Type 


Effective Four Quark Operator 




16m|^ L 

+ {dk-ff,Lui) (ua^Ldj) + 
+ 2 (dki^Ld,) (dn^Ldj) 


(dk^^Ldi) (ua^'Luj) 

(uk^^Ldi) {da^Luj) 
+ /i.e. 


$3 


-^-^ — [Uk'j^j.Lui) [da'^Ldj) + {dkj^Ldi) [ua'^Luj) 

- {dk^fj,Lui) {ua^Ldj) - (uk'j^Ldi) (da^Luj) + h.c. 


— ^^-^ — {uk^t^Rui) {da^Rdj) + h.c. 


^^M? r - 

-g^--2 — 1^ {dkRui) (uiRdj) - (ukRui) {diRdj) 

+ liidkf^'^" Rui) (uia^uRdj) - {ukO^^Rui) {dia^pRdj)} + h.c. 




h^'h^'^ 

-^-^{dk^f^Rd,) {da^Rdj) + h.c. 


Vi 


, ^ 2 {uk^fj,Lui) + {dk^fj,Ldi) 






— -^^^ {dk-i^Ldi){ua''Ruj) + /i.e. 

4my^ 



Table 2: T/ie effective four quark operator for various diquarks. The operators for ^4:, 

V2 <ind V4 mirror those for <I>i, <I>3, $7, Vi and V3 respectively, albeit with a different colour 

factor (see text). Here L{R) = 1 ^75. 



(Ic ® lc)i and (8c ® 8c)i current structures. Thus, transforming {qi^qj)[qi^' Qk) *° ^he 
form {qk^" Qi){Qi^"' Qj) implies that we are dealing with linear combinations of the form 

2 

(3c0 3c)i = -(lcOlc)i - (8c0 8c)i 

(4) 

(6c8)6c)i = -(lc^lc)i + -(8c®8c)i 

These extra color-factors need to be included while calculating the hadronic matrix ele- 
ments. 

Looking at Table it is obvious that the effective Hamiltonian for the full theory can be 
parametrized as 

9 



=0 
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with 

Ho = 

n2 = 

He = 



{qi-/f,Rb){q2-/f'Rq3) 
{qil,,Rb){q2l^Lq'i) 
{qiLh){q2Rqz) 
(giL6) (52^93) 



Hi = 
Hs = 

= 

Hj = 



(gi7^L6)(g27^^93) 

(gi7^L6)(g27^^93) 
{qiRb){q2Lqs) 

iqiRb){q2Rq3) 



(5b) 



Hs = {qia^,Lb){q2a''''Lq3) = {qia^,Rb){q2a^'^ Rq^). 

The strengths bi include both the SM contributions as well as diquark contributions (as in 
Table ^) wherever applicable. 

It now remains to calculate the hadronic matrix elements for Tii, a task that is rendered very 
difficult by the associated strong interaction dynamics. Hence, instead of attempting an exact 
calculation, one normally takes recourse to some appropriate approximation. In the "Naive 
Factorisation" approach [|l^] , the matrix elements of a four-quark operator are approximated 
by products of matrix elements of the associated quark bilinears. As an example, the amplitude 
for the decay B Xi + X2 (where Xi^2 are arbitrary mesons) can be expressed as 

{XiX2{\qirb){q2r'q3)\B) « {X^\qirb\B) {X2\q2r' q-M + {X2\qirb\B) {X^\q2r' q^lO) 

+ ^ {Xi\qiT"b\B) {X2\q2r"'q3\0) + ^ {X2\qir"b\B) {Xi\q2r'"q3\0) 

(6) 

where the second line refers to Fierz rearranged currents. Of course, only some of the matrix 
elements are non- vanishing. For one, within this approximation, the contributions of the tensor 
operators in eq.(^) vanish identically. It should be noted that eqs.(|5b| & ^) contain only color- 
singlet currents. For color-octet currents to contribute, one would need to consider additional 
gluon exchanges. Such effects are clearly not factorizable. Within this approximation then, the 
color octet parts of eq.(^) can be neglected, or, in other words, the contribution of a color-sextet 
diquark is almost indistinguishable from that of the corresponding color-triplet one. 



3 Hadronic B decays 

Within the SM, hadronic B decays may proceed through either tree level W boson exchange 
diagrams and/or through penguin diagrams (both QCD and electroweak) . The corresponding 



effective Hamiltonian, including the QCD corrections have been presented in Refs. |14,p!q|. 
For brevity's sake, we do not repeat the entire list here. It suffices to remember that the 
SM amplitudes are proportional to the Fermi constant Gp, the relevant product of two CKM 
matrix elements VibV*^ and/or VtbV^^ (with i and j as generic up type quarks and k as down 
type quark) and the combination of the Wilson coefficients that incorporates the short distance 
QCD corrections at the B mass scale. The considerable suppression due to the smallness of 
the CKM mixing is what makes -B-decays sensitive to new physics effects. 

Reverting to the calculation of the hadronic matrix elements, the decay constant fi for a 
generic (pseudoscalar or vector) meson is defined through the relations 

{P{pp)\qjipii^qi\^) = -ifppp 

{Vipv)\qjl^.qi\0) = fymve^'. ^ ' 

Here it is assumed that the meson is composed of a qjqi pair. The decay constants are best 
determined from an analysis of the respective leptonic decay modes and the relevant ones 
are listed in Table ^. The matrix elements for the associated density operators may then be 
evaluated using the Dirac equation: 

d^iqaaqj) = i{mj-mi)qiqj 
d'^iqUalbqj) = i{mj +mi)qi^5qj. 
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u 


fp 


fK 


fK' 


fo 


fn- 






fj/-f 


131 


207 


158 


214 


200 


230 


250 


275 


405 



Table 3: Values of Decay Constants in MeV. 

The matrix elements for quark bilinears between a B meson and a pseudoscalar /vector 
meson can be parametrized in terms of form factors: 



{P{pp)\qjl^{l--f5)b\B{pB)) 
{V{pv)\ni^-l5)b\BipB)) 



(PB +Pp)f, 



m 



B 



mi 



-e/xi^Q/je PbPv', 



{mB + my) 



'-qf.F^{q^ 



'-qfi){'mB + mv)Ai{q'^ 



+ i\iPB+Pv) 



(ml - M2 ) 



Qp 



A2iq' 



{tub + my) 



: 2mv{e*-g) 



q/iA-oiq'^ 



(8) 

where q = Pb — Pp(v) ^ polarisation vector of y. The apparent poles at = are 

fictitious since 

Fi(0) = Fo(0) 
2mvAo{0) = {niB + mv)Ai{0) - {niB -mv)A2iO) . 

The numerical values of the form factors can be calculated within a given model. For our 



analysis, we adopt the BSW model |16,17], and the relevant form factors, at zero momentum 
transfer, are given in Table |^ jl^,^]. It can easily be checked that choosing a different model 
for the calculation of hadronic matrix elements would not change our results appreciably. 



Decay Mode 


Fi(0) 


i^o(O) 


V{0) 


^i(O) 


^2(0) 


^o(O) 


B ^TT 


0.33 


0.33 










B ^ K 


0.38 


0.38 










B ^ D 


0.69 


0.69 










B^p 






0.33 


0.28 


0.28 


0.28 


B K* 






0.37 


0.33 


0.33 


0.32 


B ^ D* 






0.71 


0.65 


0.69 


0.62 



Table 4: Form Factors at Zero Momentum Transfer in the BSW Model 



For the q^ dependence of these form factors we assume a simple pole formula |14, 17] 
F{q'^) = F(0)/(1 — q^ / tn^joie) with the pole mass rupoie the same as that of the lowest lying 
meson with the appropriate quantum numbers (J^ = O"*" for Fq] 1~ for Fi and V; l"*" for Ai 
and A2\ 0~ for Aq). The values of these pole masses are presented in Table ^ [|l4, 17|. 
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Current 


m(O-) 


m(l ) 


m(l+) 


m(0+) 


ub 


5.28 


5.32 


5.37 


5.73 


db 


5.28 


5.32 


5.37 


5.73 


sb 


5.37 


5.41 


5.82 


5.89 


cb 


6.30 


6.34 


6.73 


6.80 



Table 5: Values of Pole Masses in GeV. 



With eqns.(0 & in place, calculation of the full matrix elements, within the factorisation 
approximation, is now a straightforward task. Consider the decay B{bq4) — > Mi{qiq4^)M2{q2q3) 
where Mj are generic mesons (pseudoscalar or vector). For simplicity's sake, assume that no 
two quarks are identical so that the quark bilinears (see eq.^) cannot relate the B to M2. In 
this case, the amplitudes are given by 



A[B- 


^PlP2] = 


A[B- 


-^PlV2] = 


A[B- 


-^ViP2] = 


A[B - 


-^ViV2] = 



[m 



P2' 



-bo + 61 + 62 



(64 - &5 - fee + fe7)"ip2 
{nib - mq^){mq,^ + m, 



2fv,mv, F^^P{ml^) [60 + 61 + 62 + 63] (e* • pp, 
'2mvjp2 (e* -pp^) 



60 + 61 - 62 - 63 + 



b^ + b^-be- b-j)m1,^ 



fv2'mv2 



fi* hi a 13 



{niB + 



(mb + mqj{mq2 + lUq^) 

(60 + 61 + 62 + 63) 



-i{el ■ €2){mB + mvi)Ai{my^) (-60 + 61 - 62 + ^3) 
+2i{pB ■ e2){pB ■ '\ . (-60 + 61 - 62 + 63) 

(9) 

For decay modes wherein qz = qi, the second set amplitudes in eq.@ contribute too. These 
additional pieces, however, can be easily read off from eq.(^. 



4 50 - 50 Mixing 

The main motivation for considering B^ — B^ mixing to constrain diquark couplings is that 
this mixing is mediated by flavour changing neutral current, which is forbidden at the tree 
level in SM. The mixing is characterised by the experimentally measurable mass difference 



i\ivid — tub^ — mB^ — 



a 



(10) 



with H and L denoting heavy and light mass eigen states. The recent world average value of 
AMd at la limit is 

AMd = (0.472 ± 0.017) x lO'^^s'^. (11) 



In SM, this mixing proceeds through the box diagrams with internal top quark and W 
boson exchanges |1^, 2^. The diagrams in which one or both top quarks are replaced by 



8 



up or charm quarks are negligible on account of: (i) the small mixing angles and (ii) the 
corresponding loop integrals being suppressed to a great extent due to the smallness of the 
light quark masses. Integrating out the internal particles, one thus gets an effective four quark 
interaction, with a {V — A)(^{V — A) current structure, and scaling as to m'y^,G'p\VtbV*f^\'^ . The 
short distance QCD corrections are well determined |T9| , pO[| , while the long-distance corrections 
are estimated to small, unlike in the case oi K — K mixing. 

In presence of diquarks, two different types of contributions may appear. If there exist 
A6 = Ad = 2 operators, then such a mixing can occur at the tree level itself. Else, additional 
contributions may appear in the form of new diquark-mediated box diagrams. However, as we 
are interested in small diquark couplings, we shall confine ourselves to tee level (in diquarks) 
processes only. 

In the calculation of the hadronic matrix element {B^\Hi\B^, the vacuum saturation ap- 
proximation is a convenient one. Herein, one inserts a complete set of states between the 
two currents and assumes that the sum is dominated by the vacuum and thus the hadronic 
matrix elements are proportional to by virtue of eq.(0). The bag factor, Bb, introduced 
to parametrize all possible deviations from the vacuum saturation approximation, can be eval- 
uated in various nonperturbative approaches. We use here the values of BsifJ-b) and fs as 
obtained by UKQCD collaboration in a quenched lattice calculation ||2T|| . 

Incorporating the contributions of all such current structures, in addition to that of the 
SM, we obtain 



AMa = flBBTUB, 



bQ + bi-b2-h + (^4 + h 

(nib +mdY 



(12) 



This may then be compared with the experimental value to obtain the required constraints. 



5 Results 

Before we determine the bounds obtainable from i?-phenomenology, it is worthwhile to reex- 
amine the SM predictions for the decay modes of interest; this helps in selecting the channels 
likely to result in stronger constraints. As mentioned earlier, within the SM, the hadronic 
B decays are mediated by one or more of tree (VF-mediated), electroweak penguin and QCD 
penguin diagrams. The branching fractions are determined primarily by the CKM mixings 
operative in the particular decay and, in case of one-loop processes, by the corresponding loop 
integral. For example, the decays B~ — > D~tt^ and B^ — > -K~7r~^ are suppressed in compar- 
ison to B'^ — > D~^7T~ , D~ , B~ — > D^D~ and B'^ — > -D^vr" on account of Vub being much 
smaller than Vcb- Similar statements obviously hold for the decays into the corresponding 
excited states. Inspite of such a suppression, the tree diagram far outweighs the one-loop 
contributions for any of these decays. For the decays B~ K^ir^, B^ K^i:^ (and the 
corresponding PV and VV modes) though, the tree level contributions are double Cabibbo 
suppressed with the consequence that these decays are dominated by penguin diagrams 0. De- 
cays like B~ — > K~K^, B~ tt'K^, on the other hand, are governed solely by electroweak 
and/or QCD penguins. 

A different suppression occurs for the processes B^ — > vr'^vr'^, B^ — > D^tt^. Compared to 
the analogous modes B^ — > 7r"'~7r~, B^ D~^tt~ wherein the charged mesons are created 
from the vacuum by a color-singlet current, these decays are obviously color-suppressed. In 



fact, the respective short distance coefficients differ by as much as a factor of 20 [14|. And 
finally, there are the annihilation diagrams in the decays like B{bq) — > Xi{qiq)X2{q2Qi)- In 
these decays, b and q in B meson annihilate to produce q and q2 quarks, which, in turn 



In our numerical calculations, we have used the Wilson coefficients as listed in Ref. |14| for Nc = 3. 
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form final state mesons with qiqi pair, created from vacuum. As these contributions are 
proportional to the wavefunction at zero, they are typically much smaller than either of the 
tree or penguin mediated spectator contributions. For example, Ref. |14] argues that the 
annihilation amplitude for B — > P1P2 modes is proportional to the mass difference of the 
mesons in the final state and hence there is essentially no annihilation contribution to S — > 
vrOvrO, K+K~ etc. 

It is tempting to assume that the modes suppressed within the SM would be the ones 
most sensitive to effects from new physics. While this is largely so, a few points should 
be remembered. For one, color suppression and/or suppression of annihilation diagrams are 
essentially independent of weak matrix elements, and hence equally applicable to either the SM 
or a theory with diquarks. Secondly, even for decays wherein the SM amplitudes are Cabibbo 
suppressed, the experimental data may not be precise enough for it to be a very sensitive 
probe. Rather, it could well turn out that an unsuppressed mode may turn out to be one of 
the most sensitive on account of the observations matching very well with the SM predictions. 

In obtaining numerical results, we assume that only one pair of diquark couplings are 
nonzero. While this restriction may seem unwarranted, it is an useful approximation that 
allows one a quantitative appreciation of the various experimental constraints. Furthermore, 
we assume a common mass of lOOGeV for all the diquarks. As explained earlier, mass splittings 
between states in a single multiplet is disfavoured by LEP data. And since the effective four- 
Fermi operator goes as rn'^'^y^, for a general diquark mass, all our bounds on the products need 
only be rescaled by a factor of (m$(y)/100 GeV)^. 

The bounds, as obtained from a given decay mode, can be broadly classified into two sets. 
An experimentally observed channel (with an associated error bar) would, in general, allow 
the diquark coupling pairs to lie in one of two non-contiguous windows, with the separation 
between the windows determined by the agreement of the SM contribution. On the other hand, 
decay modes that are yet to be experimentally seen, can only lead to a single window. For 
a specific combination of diquark couplings, we look at all such individual bounds and then 
delineate the range satisfied by each. As an illustration, let us consider the product h^i-]h^i2 ■ 
At 90% C.L., the ranges allowed by individual decays are as follows: 



(a): 


B'> - 


D+p- 


[-6.0 X 10" 


-2 


-5.3 X 10-2], 


[-7.3 X IQ-'^ 


5.7 X 10"^ 


(6): 


B- - 


D'^p- 


[-4.7 X 10" 


-2 


-4.2 X 10-'^], 


[-4.7 X lO"'' 


4.8 X IQ-^ 


(c): 


B- - 


K°*7r- 


[-8.9 X 10" 


-4 


2.2 X IQ-^] 






id): 


B^ - 


K^p'' 


[-2.1 X 10" 


-3 


1.6 X lO^^] 







Clearly the first window allowed by (a) is completely ruled out by the the bounds from (c) 
and (d). The same is true for the first window from (6). Progressively eliminating parts of 
the domains allowed by the individual decays, we find that the actual allowed range for this 
particular combination is only [—4.7 x 10~'^,5.7 x 10^^]. An identical strategy is adopted for 
all other combinations, and we list the best bounds in Tables |6HT3|. 

A very important point is to be noted here. In the preceding analysis, while we have 
selected the range of parameters common to each constraint, we have not really used the entire 
information available to us. Such an analysis would involve the use of a statistical discriminator 
such as a test or a likelihood test. While such an exercise is a straightforward one and would 
have led to bounds stricter than those we list, the decision to forego it was a conscious one. 
For, in the absence of higher order corrections and a more precise calculation of the hadronic 
matrix elements, the bounds derived here are only indicative. Hence, further refinement using 
statistical methods is not really called for. 

From Table I, it is easy to see that /.«4^\ /^Jf /ig^, h^h^, /.Jf h^h^J^^ and /^Jf 4? 
result in both neutral and charged current structures. In a given hadronic decay, both the op- 
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erators contribute with one of them being color suppressed. Furthermore, the flavor structure 
determines whether two contributions interfere constructively or destructively. As a conse- 
quence, the bounds on hf^h^^^ , ^if^ld ^^'^ weaker by a factor of compared 
to those for h^tj^h^^^ and h^ij^h!~^i . A color-unsuppressed operator associated with com- 
binations h^^j^h!'^^ and h^tj^h^^^ are neutral current ones, their contributions to charged current 
decays are naturally color-suppressed. Hence the corresponding bounds are weaker. Similarly, 
since h[^^h^^j^ and h^'^^h^^j^ , are associated only with scalar, pseudoscalar and tensor operators, 
they cannot contribute to i? — > VV decays. And finally, as the diquarks ^7^8 couple only 
to down-type quarks, they can contribute only to those decays that occur in the SM solely 
through penguin diagrams. 

The diquark $3 differs from only by colour quantum number. But since the color 
factors for triplet and sextet diquarks are accidentally equal and h^^^ and h['^^ 's have no specific 
symmetry property under the exchange of i and j, the bounds on the product of these couplings 
are exactly the same. A similar story obtains for other sets of diquarks, (Vi, V2) and (^3, V4). 

As discussed earlier, h\j^ is analogous to the trilinear R parity violating coupling X'-jf,. Thus 
the constraints on h[j^h^j^j^ are equivalent to those on A-'^^A^'^;. The upper bound on h!2z'^^22 
is marginally weaker than that, quoted in |12| and |22|. The upper bound on 'tiiifi'i^ is much 
weaker than that listed in [^] and ||2^ . For other h^ij^li'^^, we obtain more stringent bounds. 
For most of the nonsupersymmetric product of two diquark couplings, our predicted bounds 
are much stronger. 

In a few hadronic decays {B~ — > K~ J/^, B~ tt~J/^ etc.), the SM predictions do not 
agree with the experimental observation even at 2a level. While this could be construed as an 
indication of new physics (in our case diquarks), we prefer to tread a more conservative path. 
Consequently, we have not included such decays in our analysis. 

Turning to B^-B^ mixing, it is obvious that a tree-level contribution will accrue from any 
four-quark operator that violates both b- and d-number by two units each. Furthermore, ^2, 
<I>3 and ^8 do not contribute to B^ — B^ mixing by virtue of the antisymmetric nature of their 
couplings under exchange of the flavor indices. Thus B'^ — B^ mixing only imposes limits on 
the parameter space of ^>i, $7, Vi and V2 diquarks. Unhke in the SM, B^ - B^ occurs at tree 
level in presence of nonzero diquark couplings. Accordingly, we obtain most stringent bounds 
on hifh[f and ^if'&if. 



6 Conclusion 

In summary, we have studied the leading order effects of scalar and vector diquark and/or R 
parity violating couplings on hadronic B decays and B—B mixing. Clearly a diquark must have 
more than one non-zero couplings to SM fields to be able to mediate such processes. We take the 
economic standpoint that only any two of such couplings are non-zero. Analysing the present 
data on 5-decays, we derive constraints on such pairs that are significantly stronger than those 
derived from other low energy processes. Theoretical improvements on nonfactorisation effects 
and estimates of annihilation form factors as well as precise measurements of the decay modes 
at the upcoming B factories in near future will improve our bounds on the parameter space 
for diquarks and/or R parity violating couplings in the minimal supersymmetric SM. 
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Product of Couplings 


Mode 


Allowed Region 






[-2.6 X 10-^l.l X 10-3] 


lo Iz 


B- D^p-, B K^p^ 


[-4.7 X 10-^,5.7 X 10^4] 


''23 "'22 


BO D+D-, 5° ^ D+*D- 


[-5.5 X 10-2,-5.3 X 10-2], 

r Q \/ 1 n~3 o o w 1 n~3i 
[ — D.o X lU , Z.Z X iU J 


^23^ ^\.2 


BO ^ i^o^o 


[-1.4 X 10-^9.5 X 10-^] 


"33 "11 


Bj - Bj Mixing 


[-1.4 X 10-^,-1.3 X 10-^], 
[-1.3 X 10-^2.4 X 10-9] 



Table 6: Bounds on $i couplings in units of (m$^/100 GeV^ at 90% C.L. 



Product of Couplings 


Mode 


Allowed Region 


"13 "l2 


B- K'TT^, B- ir-K^, 
B- K-p^ 


[-1.2 X 10-2,7.7 X 10-"^] 


^23^ ^12^ 


BO ^ i^o^o 


[-2.8 X 10-^ 1.9 X 10-3] 



Table 7: Bounds on $2 couplings in units of fm^j/lOO GeVj at 90% C.L. 
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Product of Couplings 


Mode 


Allowed Region 


"is "'12 


B- ^TT^K-, B- ^K-p^, 


[-1.9 X 10~^,2.9 X 10"^] 


"23 12 


B" -s- 7r°Dj, 5° ^ /9°J/* 


[-1.5 X 10~^, 1.7 X 10~^] 


"13 "11 




[-3.3 X 10~^, 2.9 X 10"^] 


'''13 "'12 




[-1.1 X 10-^7.8 X 10"^] 


^(3)^(3) 
"'13 "21 




[-1.2 X 10~^, 1.4 X 10~^] 


"23 "22 


B- D^D-*, 5° ^ D+D- 


[-8.7 X 10-3,3.1 X 10-2] 


^^23^ ^21^ 




[-4.2 X 10-3,4.2 X 10-3] 


^23 ^12 




[-1.4 X 10-2,2.0 X 10-2] 


13 11 


50 ^ TtOttO, 5° ^ TT+TT- 


[-5.1 X 10-3^7 4 X 10-3] 


"l2 "l3 




f— 1 4 X in~3 9 9 X in-3i 


"l3 "l2 


B- ^K'tt", B^^tt+K-, 


[-3.3 X 10-3, 1.4 X 10-3] , 


"l3 "21 


B- D^p- 


[-5.2 X 10-3,2.2 X 10-3] 


"23 "'22 


BO ^ D+DJ, B- ^ Z^OZ),-, 


[-4.9 X 10-3,1.5 X 10-2] 


^12'' ^23'' 




[-1.6 X 10-2,1.1 X 10-2] 


^23^ ^12^ 


B- ^ TT^i^,- 


[-1.1 X 10-2,1.6 X 10-2] 



Table 8: Bounds on $3 couplings in units of GeV) at 90% C.L. 
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Product of Couplings 


Mode 


Allowed Region 


""IS ''11 




[-1.1 X 10-3,2.6 X 10-2] 


"13 "'12 




[-8.1 X 10-^,5.4 X lO-''] 


"23 "22 


^ ^ D+*D- 

B- D^D- 


[-2.2 X 10-3,6.5 X 10-3], 
[5.3 X 10-2,5.5 X 10-2] 


^23 ^12 


B- 7r°Dj,50 ^ pOj/* 


[-3.4 X 10-3,3.1 X 10-3] 


"l3 "ll 


fiO ^ 7r°7r°, ^ tt-ttO 


[-3.3 X 10-3,2.9 X 10-3] 


^(4);^(4) 
"l3 "l2 




[-1.1 X 10-3,7.8 X 10-^] 


K(4)r(4) 

"13 "21 


B- TT-D^, B- p-D^, 


[-1.2 X 10-3, 1.4 X 10-3] 


K(4)^(4) 
"'23 22 


B~ D^D~*, B° D+D~ 


[-8.7 X 10-3,3.1 X 10-2] 


I(4)r(4) 

"23 "21 




[-4.2 X 10-3,4.2 X 10-3] 


^12^ 




[-1.4 X 10-2,2.0 X 10-2] 


11 13 


B^^TT+TT-, SO^TtOttO 


[-3.6 X 10-3,6.4 X 10-3] 


"13 "ll 




[-6.4 X 10-3,3.6 X 10-3] 


"l2 "l3 


B^ K+K-, B° ^ TT+K- 


[-1.5 X 10-2,-1.2 X 10-2 
[-2.3 X 10-3,1.5 X 10-3] 


"l3 "l2 


^0 ^ ^0*^0^ B- K-p^ 


[-1.5 X 10-3,2.3 X 10-3}] 


"13 21 


BO ^ ^0^0^ B- ^ L>V- 


[-3.7 X 10-3,2.2 X 10-2] 


"22 "23 




[-3.5 X 10-3,1.0 X 10-2] 
f7 4 X 10-2 w -|n-2i 


"23 "22 


B^-.D+D-,B- ^D^D-, 
BO D+*D- 


[-1.0 X 10-2,3.5 X 10-3] 




B- ^ -"17- 


[-8.0 X 10-3.1.2 X 10-2] 


^23^^12^ 




[-1.2 X 10-2,8.0 X 10-3] 



Table 9: Bounds on $4 couplings in units of f'm$4/100 GeVj at 90% C.L. 
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Product of Couplings 


Mode 


Allowed Region 


"■13 'hi 




[-1.1 X 10-3,1.3 X 10-3] 


'ns "'12 


B- TT-R^, B- TT^K- 


[-3.9 X 10-^,7.3 X 10-^] 


^23 ^12 




[-9.5 X 10-^, 1.4 X 10-3] 


"33 "ll 


- Mixing 


[-1.4 X 10-^-1.3 X lO"'^], 
[-1.3 X 10-^2.4 X lO-'^] 



Table 10: Bounds on $7 couplings in units of GeV^ at 90% C.L. 



Product of Couplings 


Mode 


Allowed Region 


"13 "12 


B- ir-K^, B- TT^K- 
B- 7r-^°* 


[-7.9 X 10-^, 1.2 X 10-3] 


^23^ ^12^ 


50 ^ ^0^0 


[-1.9 X 10-3,2.8 X 10-3] 



Table 11: Bounds on couplings in units of (m$g/100 GeVj at 90% C.L. 
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Product of Couplings 


Mode 


Allowed Region 


^(1)^(1) 

"13 "11 


BO ^ ;r+7r-, ^ pOp° 


[-6.2 X 10-3,3.2 X 10-3] 


'^31 '^11 




[-2.1 X 10^3^1.7 X 10~3] 


^13 ^12 




[-1.4 X 10-3,2.0 X 10^3]^ 
[1.0 X 10-2,1.1 X 10-2] 


'^13 "21 


B- D^p- 


[-1.0 X 10-3,-2.6 X 10-3], 
[-8.1 X 10-'',9.8 X 10-^] 


^^31 ^12 


BO ifOyrO, B- ^ TT-^o 


[-9.8 X 10-'',8.1 X 10-^] 


^^31 "21 


^ 7r°if-, B- p^K- 


[-1.5 X 10-3,5.3 X 10-^] 


^^23 ^^22 


BO _^ D+D-, B- D^DJ 
B- D^*D- 


[-9.2 X 10-3,3.1 X 10-3] 


^^23 ^12 


BO _ i^o^o 


[-1.3 X 10-3,1.9 X 10-3] 






[-2.1 X 10-3.2.1 X 10-3] 


^^32 "12 


BO ^ ^0^0 


[-4.7 X 10-3,7.0 X 10-3] 


^^32 '^21 




[-1.9 X 10-3,1.3 X 10-3] 


^^33 '^ll 


5^ - Mixing 


[-5.4 X 10-^,-4.8 X 10-^], 
[-4.9 X 10-9,8.9 X 10-^°] 



Table 12: Bounds on V\ couplings in units of (mvi/WO GeV) at 90% C.L. 
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Product of Couplings 


Mode 


Allowed Region 


^(3).{3) 




[-3.2 X 10-^1.0 X 10-2] 


^(3) .(3) 
''31 ''11 




[-1.4 X 10-3,1.8 X 10^3] 


^(3) .(3) 
'^13 '^12 


7T+K-, 5° ^ p+K- 


[-1.2 X 10-2,-1.0 X 10-2], 
[-2.0 X 10-2, 1.4 X 10-'^] 


"13 "21 




f-1 X 10"^ 9 X IQ-^l 


^(3) .(3) 
^31 ^12 




[-1.5 X 10-^4.1 X 10-3] 


^(3) .(3) 
''31 '^21 


B- ir^K-, B- K-p° 


[-7.3 X 10-^1.7 X 10-3] 


^(3)^(3) 
''32 ''22 


BO D+D-, B- D^D~ 


[-1.5 X 10-3,4.3 X 10-3] 


^(3)„(3) 
'^23 '^22 


B- D^*D- 


[-2.2 X 10-\-1.4 X 10-1], 
[-2.0 X 10-2,5.7 X 10-2] 


^(3)^(3) 
'^23 '^12 




[-1.1 X 10-2,1.4 X 10-2] 


^(3)^(3) 
'^23 '^21 


B^ p^J/^ 


[-2.5 X 10-3,2.1 X 10-3] 


^(3)^(3) 
''32 ''12 


5° ^ TfO J/* 


[-2.5 X 10-3,2.1 X 10-3] 


^(3). (3) 
'^32 '^21 


5- Z^.-TfO 


[-7.0 X 10-3,1.0 X 10-2] 



Table 13: Bounds on V3 couplings in units of (myg/lOO GeV^ at 90% C.L. 
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